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Abstract 

The present matter density of the Universe, while highly inhomogeneous on small scales, 
displays approximate homogeneity on large scales. We propose that whereas it is justified to 
use the Friedmann-Lemaitre-Robertson- Walker (FLRW) line element (which describes an 
exactly homogeneous and isotropic universe) as a template to construct luminosity distances 
in order to compare observations with theory, the evolution of the scale factor in such a 
construction must be governed not by the standard Einstein equations for the FLRW met- 
ric, but by the modified Friedmann equations derived by Buchert [H [8] in the context of 
spatial averaging in Cosmology. Furthermore, we argue that this scale factor, defined in the 
spatially averaged cosmology, will correspond to the effective FLRW metric provided the size 
of the averaging domain coincides with the scale at which cosmological homogeneity arises. 
This allows us, in principle, to compare predictions of a spatially averaged cosmology with 
observations, in the standard manner, for instance by computing the luminosity distance 
versus red-shift relation. The predictions of the spatially averaged cosmology would in gen- 
eral differ from standard FLRW cosmology, because the scale-factor now obeys the modified 
FLRW equations. This could help determine, by comparing with observations, whether or 
not cosmological inhomogeneities are an alternative explanation for the observed cosmic 
acceleration. 



1 Introduction 

The goal of Cosmology is to describe the evolution history of the Universe as a whole, tracing it 
as far back as possible from the present epoch. This seemingly daunting task is made tractable 
by the fact that the observed matter density in the Universe appears to be approximately 
homogeneous when viewed on the largest scales. Further, observations of the Cosmic Microwave 
Background (CMB) radiation reveal it to be isotropic (after subtracting a dipole contribution 
due to our local motion) down to the level of one part in 10^ . The twin conditions of homogeneity 
and isotropy, if valid on all length scales, would immensely simplify the form of the metric and 
matter tensors used in the General Relativistic description of Cosmology. It can be shown 
[T] that on purely geometric grounds, the most general metric describing a homogeneous and 
isotropic universe takes on the form of the Friedmann-Lemaitre-Robertson- Walker (FLRW) line 
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element. The standard method of proceeding is to assume that the geometry of the universe 
can be described by a background metric with small perturbations, where the background is 
exactly homogeneous and isotropic and is hence of the FLRW form. (For a detailed discussion 
of the assumptions involved in the standard Cosmology, see Ref. [2].) 

This approach is adopted not only for the early universe, for example in describing Big 
Bang nucleosynthesis and the dynamics of the CMB anisotropics, but is also assumed to be 
valid at recent epochs. In particular, the luminosity distances constructed to compare data from 
Type la Supernovae (SNe) with theory, assume a metric of the FLRW form with its dynamics 
given by the Einstein equations for a homogeneous and isotropic universe. We know, however, 
that the matter content of the Universe today is not homogeneous on small scales - one finds 
inhomogeneity at the scale of individual galaxies for example; and there is evidence |3j that 
voids of order 30h' in diameter account for 40-50% of the observed Universe. There is 

also evidence for voids 3-5 times this size [4j, as well as local voids on smaller scales [5]. The 
evidence for homogeneity of the matter distribution in the Universe only appears at scales larger 
than about lOO/i"^ Mpc or so p] . It is then reasonable to assume that the metric describing our 
Universe must also reflect the inhomogeneity of the matter on the corresponding scales. One 
may now ask whether it is at all justified to use the luminosity distances constructed using the 
FLRW metric with its associated dynamics to compare theory with observational data from SNe. 
In this paper we argue that the use of such luminosity distances is indeed justified, provided 
we are only interested in using such constructions in probing length scales larger than a certain 
minimum scale which is fixed by requiring that the assumption of homogeneity be valid (a 
lower bound on such a scale could be 100/i~^ Mpc for example). More importantly, we propose 
that the dynamics of this (effective or template) macroscopic metric must be governed not by 
the standard Einstein equations applied to the FLRW metric, but by the effective equations 
obtained after spatial averaging on constant time slices, as developed by Buchert 018]. 

An outstanding open issue has been the following: how does one compare the predictions of 
the inhomogeneous, spatially averaged cosmology with observations, and should one be attempt- 
ing an interpretation in terms of spatially averaged variables in the first place? For instance, one 
could study light propagation and determine the luminosity distance versus redshift {Di^{z)) 
relation which best fits the observations, within an inhomogeneous model of the Universe. Fol- 
lowing the important and interesting early work of Dyer and Roeder, where the angular-size 
and luminosity distances were obtained for beams propagating in "lumpy" FLRW universes 
[9], detailed investigations have been carried out in [lOj. Such a program has also received 
considerable attention in the context of the spherically symmetric inhomogeneous Lemaitre- 
Tolman-Bondi (LTB) models [11^ I12j . However, it is not clear to us that this is necessarily 
the best way to proceed, at least as far as distant SNe are concerned. A realistic computation 
of light propagation in an inhomogeneous geometry seems difficult at best, and while the use 
of simplified toy models such as the LTB solutions in this endeavor may be justified on the 
grounds that we may inhabit a locally underdense region [4], it is not clear that the assump- 
tions justifying the use of the toy model remain valid at large scales (see also Ref. |12j for 
a recent review of the LTB approach). And after all, the concordance model infers the pres- 
ence of a dark energy by comparing the predictions of the standard homogeneous and isotropic 
FLRW cosmology against observations [T3]. If one decides to interpret the FLRW cosmology 
as a template against which observations may be compared, then in order to decide whether 

^Here h is the conventional parameter which appears in the definition of the Hubble constant Ho — WOh 
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cosmic inhomogeneities could be an alternative to the concordance model, one must construct 
a homogeneous and isotropic cosmology, after a suitable smoothing procedure, and compare its 
predictions with observations such as the Di^{z) relation for Type la SNe. 

One should note that the effects of averaging within the standard perturbed FLRW frame- 
work of Cosmology, have been studied by several authors (see, e.g., Refs. |14l I15| ). and it has 
been argued that these effects must necessarily be small. The papers in Refs. [14^ I15j however, 
do not necessarily address the entire problem. For example, Ref. [Hj deals with only super- 
horizon fluctuations, whereas work done by others (see, e.g, Ref. [16]) indicates that there may 
be nontrivial contributions due to averaging from subhoiizon fluctuations where the matter 
perturbations have gone nonlinear. And while the authors of Ref. [15] construct a cosmo- 
logical model with negligible effects of averaging, there do exist toy models of inhomogeneous 
spacetimes in which these effects are large (see, e.g., Ref. [17\). 

Furthermore, the possibility of nonperturbative effects of averaging has not, to our knowl- 
edge, been satisfactorily addressed, much less ruled out. It is important to realize that current 
A^-body simulations for example, which are based on a Newtonian approximation of gravity, 
consistently miss possible relativistic effects of averaging. This is because the modifications due 
to averaging in a Newtonian setting can be shown to reduce to boundary terms in general |18j . 
and to vanish in the case of spherical symmetry [19] . These terms therefore do not contribute 
in simulations which either set periodic boundary conditions on the simulation box, or use a 
spherically symmetric simulation region. As we will see, the "backreaction" effects in Buchert's 
framework (i.e., the modifications to the Cosmological equations due to averaging) are coupled 
to the average of the spatial curvature, this being a completely relativistic effect. Nontrivial 
consequences of this coupling deserve more attention than they have received in the context of 
A^-body simulations (see Ref. [T9]). 

We also note that an argument for describing the Universe using FLRW plus Newtonian 
perturbations (e.g., as given in Ref. [2D]), is based on the fact that it is Newtonian perturbations 
and a negative pressure component, which describe the Universe very well. In other words, one 
argues that the effect of averaging inhomogeneities must be small, but at the cost of adding 
an ill-understood Dark Energy component in the Universe. On the other hand one should 
allow for the following possibility (as of now neither confirmed nor ruled out): at an initial 
epoch, say at the end of inflation, the perturbations in the homogeneous metric would give a 
small "backreaction" as a result of averaging. As a result of its coupling with the curvature, 
this backreaction would grow and become signiflcant around the time of structure formation, 
when matter inhomogeneities on small scales become nonlinear. The backreaction would also 
nontrivially affect the growth of the metric perturbations. Such a possibility could only be 
tested within a fully nonperturbative treatment of the evolution of inhomogeneities. 

The goal of this paper is modest. We argue that the scale factor, deflned in the spatially 
averaged cosmology via the volume of a domain, should correspond to an effective FLRW 
metric, provided the size of the averaging domain coincides with the scale at which cosmological 
homogeneity arises. This allows us, in principle, to compare predictions of a spatially averaged 
cosmology with observations, in the standard manner, for instance by computing the Di^{z) 
relation. (However, see Ref. |21j for an alternative approach to matching the spatially averaged 
dynamics with observations.) The predictions of the spatially averaged cosmology would in 
general differ from standard FLRW cosmology, because the scale-factor now obeys modified 
Einstein equations. We will demonstrate that consistent models of spatially averaged FLRW 
templates exist which can, at least in principle, explain observations such as those of the Di^{z) 
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relation for Supernovae. 

In a nutshell, our new results in this paper are that we propose a template metric which 
takes into account the kinematical effects of inhomogeneities in the Universe. To implement 
this idea, we first suggest that in Buchert's spatial averaging scheme, the size of the averaging 
domain should coincide with the cosmological scale at which homogeneity sets in. We then argue 
that, as a first step towards predicting the relation between inhomogeneities and observations, 
the expansion factor ax)(t) defined by Buchert should be identified with the scale-factor a{t) 
in standard FLRW cosmology. The template metric is then the FLRW metric, with ax>{t) 
playing the role of a{t). We then show, with the help of a few toy examples, that the luminosity 
distance - redshift relation in a model Universe with an accelerating ax)(t) can be computed 
in the standard manner. To the best of our knowledge, such a template metric has not been 
proposed in earlier work by other researchers. 

2 The Coarse Grained Picture and Cosmology 

At the length scales of cosmological interest, individual galaxies are small enough to be treated 
as point-like objects. Given the large number (~ 10^^) of galaxies estimated to exist in the 
observed Universe, it is also assumed that the matter content of the Universe can be described to 
a good approximation as a continuum fluid. At the level of the matter distribution, this picture 
is comparable to that of a gas composed of a large number of molecules, wherein the internal 
structure of the molecules can be ignored when dealing with the macroscopic properties of the 
gas which is treated as a continuum fluid. Further, the complex intermolecular interactions in 
the gas at the microscopic level can be dealt with in a mean field approximation by considering 
volume averages of the quantities of interest. 

This picture is made precise by introducing the concept of physically infinitesimal volumes 
|22j or coarse- graining cells. The fluid approximation is assumed to be valid provided the 
(effectively) infinitesimal volumes used in defining continuum quantities such as the matter 
density of the fluid, correspond to physical volumes which are large enough to contain a very 
large number of the particles constituting the fluid. At the same time, since these volumes are 
to be effectively infinitesimal, their physical size must also be much smaller than the scale of 
the fluid as a whole. We can then say that the length scales L being probed by these physically 
infinitesimal volumes must lie in a certain rangcl Li < L < L2 which is determined by the 
details of the fluid (see also [23]). Clearly we must satisfy L2 <C Lsyst where Lgyst is the size of 
the system. In the cosmological context, Lgyst would be the Hubble scale, and in order to deal 
with homogeneous and isotropic effective descriptions we would also have to satisfy Li > L/iom 
where the homogeneity scale -L^om as mentioned earlier can be as large as 100/i~^ Mpc. The 
physical quantities we wish to deal with, such as the matter density, pressure, fluid velocity 
field, etc. are defined by suitably averaging or coarse graining over the physically infinitesimal 
volumes (hence the terminology of coarse-graining cells (CGCs)). For example, the matter 
density of a given CGC is defined as the mass contained in it divided by its physical volume. 
The (effectively) continuum matter density is built up by repeating this for every CGC in the 
fluid. 

Since we deal with the notion of a continuum fluid in Cosmology, clearly we need a working 
description of coarse graining appropriate in this context. However, the underlying theory 

^In the cosmological context, these length scales will be probed by the metric. Since we will soon specialize 
to a metric written in comoving coordinates, these length scales must be understood to be comoving scales. 
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Figure 1: Symbolic depiction of the modifications to Einstein's equations caused by an explicit 
smoothing. Adapted from Ellis [24j . fig. 7. See text for details. 



being General Relativity, there are problems. The crucial issue is that along with the matter 
distribution (which is now tensorial), the metric for the spacetime also needs to be coarse 
grained. As discussed by Ellis [23] , in defining a smoothing operation to go from the scale of stars 
(Scale 1 in Ellis' terminology, microscopic in ours) to cosmological scales where homogeneity 
can be assumed to have set in (Scale 5 for Ellis, macroscopic for us), one would employ two 
logically different smoothing operators, one each on the metric and the matter tensor. We 
denote the smoothing operator acting on the metric as Smet and the one acting on the matter 
tensor as Smatt (our notation differs from that used in Ref. [21] )■ We assume that Einstein's 
equations are valid on the microscopic scale of stars, which is reasonable since it is at these 
scales that Solar System tests and binary pulsar studies probe General Relativity. Symbolically 
then (see fig. the metric gmic at this scale is used to construct the Einstein tensor Gmic 
which is equated to the matter tensor Tmic-, giving Einstein's equations Gmic = Tmic- To obtain 
the corresponding description at the macroscopic scale, we use the smoothing operators to get a 
new metric tensor gmac = Smet {gmic) and a matter tensor Tmac = Smatt (Tmic)- The macroscopic 
metric gmac is used to construct the Einstein tensor Gmac at the macroscopic scale. Formally, 
one can think of the map Smet as inducing a map SEinstein which acts on the microscopic 
Einstein tensor Gmic to yield Gmac = SEinstein{Gmic)- Since the operators SEinstein and Smatt 
are in general different, we have Gmac 7^ Tmac- Hence, if the Einstein equations are valid on a 
given length scale, then they are in general not valid on some other length scale connected to the 
first by appropriately defined smoothing operators. If we can re-express the new Einstein tensor 
Gmac ill terms of quantities obtained by acting the operator Smatt on both sides of the Einstein 
equations at the microscopic scale, then it should be possible to write Gmac = Tmac + Pmac 
where Pmac represents corrections to the Einstein equations arising from an explicit smoothing. 
In the standard approach, it is assumed that these corrections are small enough to be neglected. 
(It should be emphasized that given the recent evidence for voids in the Universe, one or more 
additional scales may eventually have to be included in the averaging scheme.) 

Added to all this is the complication that smoothing or averaging operations are notoriously 
difficult to define for tensor quantities in a gauge covariant manner. Some examples of averaging 
techniques employed in General Relativity (and in Cosmology in particular) can be found in 
Ref. [25]. Here we restrict ourselves to the formalism developed by Buchert [3 [8] for volume 
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averages of scalar quantities in a given choice of time-slicing. In the following sections we recall 
this formalism and argue that it is naturally adapted to the notion of the smoothing operator 
for the matter sector. We will reserve the term "averaging" for the matter smoothing defined 
in this way. As for the smoothing of the metric (which we will refer to as "coarse graining" ) , we 
will take recourse to cosmological observations and argue that the result of such an operation 
on the inhomogeneous metric must yield the FLRW metric on a suitable coarse graining scale. 



3 Averaging in a Dust Dominated Spacetime 

For simplicity we restrict our attention to a pressure-less fluid or "dust", although a general- 
ization to fluids with non-zero pressure is also possible (see Ref. p]). For a general spacetime 
containing irrotational dust, the metric can be written in synchronous and comoving gaug^, 

ds'^ = -dt^ + hij{x, t)dx'dx^ . (1) 

The expansion tensor 0^- is given by 0*- = {l/2)h^^hkj where the dot refers to a derivative with 
respect to (proper) time t. The traceless symmetric shear tensor is defined as a*- = ©j — (0/3)(5* 
where = G* is the expansion scalar. The Einstein equations can be split [7j into a set of 
scalar equations and a set of vector and traceless tensor equations. The scalar equations are 
the Hamiltonian constraint (|2ap and the evolution equation for Q (I2b|l . 



(=^)7^ + - 2a^ = IGnGp , (2a) 

^^^n + e + Q^ = i27tGp, (2b) 

where ^^^TZ is the Ricci scalar of the 3-dimensional hypersurface of constant t, is the rate 
of shear defined by = (1/2)(T*(T^ and p is the matter density of the dust. Note that all 
quantities in Eqns. ([2]) generically depend on both position x and time t. Eqns. (pa|) and ()2bp 
can be combined to give Raychaudhuri's equation 

e + \e^ + 2a^ + 4ttGp = 0. (3) 

The continuity equation p = —Qp which gives the evolution of p, is consistent with Eqns. (j2ap . 
()2bp . We only consider the scalar equations, since the spatial average of a scalar quantity can be 
defined in a gauge covariant manner within a given foliation of space-time. For the space-time 
described by ([1]), the spatial average of a scalar ^'(x, t) over a comoving domain V at time t is 
defined by 

W© = T^ f d^'xVh^, (4) 
vv Jv 

where h is the determinant of the 3-metric hij and Vx> is the volume of the comoving domain 
given hy Vx> = Jj) d^xVh. Spatial averaging is, by definition, not generally covariant. Thus 
the choice of foliation is relevant, and should be motivated on physical grounds. In the context 
of cosmology, averaging over freely-falling observers is a natural choice, especially when one 
intends to compare the results with standard FLRW cosmology. Following the definition Q 
the following commutation relation then holds [7] 

- {^)v = {'^Q)v - mv{Q)v . (5) 
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which yields for the expansion scalar O 

{eyT,-{@)T> = {e')v-{e)l. (6) 

Introducing the dimensionless scale factor ax> = {Vx> /Vx)i)^^^ normalized by the volume of the 
domain V at some initial time ti, we can average the scalar Einstein equations (f2al) . (I2b|) and 
the continuity equation to obtain [7j 

{p)v= -{Q)v{p)v ; (e)i, = 3^, (7a) 

av 

3 (^)' - + \{n)v = , (7b) 

3(^)+A7rG{p)v = Qv. (7c) 



Here {TZ)xi, the average of the spatial Ricci scalar ^^^7^, is a domain dependent spatial constant. 
The 'kinematical backreaction' Qd is given by 

= ^ ((e^)^ - (6)2,) - 2{a')v (8) 

and is also a spatial constant over the domain D. A necessary condition for ([7c|) to integrate 
to (|7bp takes the form of the following differential equation involving Qd and {TV)x>, 

Qv + Q—QT> + {nyj, + 2^{n)v = o. (9) 

av av 



4 Smoothing Operators and Modified Field Equations 

The averager defined in Eqn. Q is a natural candidate for the smoothing operator Smat 
which must act on the matter tensoi0, although we have not yet specified the domain size. 
Note that this averaging is defined on an inhomogeneous underlying geometry, i.e. it does not 
smoothen out the geometry itself. It remains to determine the operator Smet that will coarse 
grain the metric. Although we do not have an explicit definition for this operator, we can 
define a notional smoothing operator acting on the metric, by using physical inputs from the 
Universe we observe. Namely, we propose that since the Universe has a matter content which 
is homogeneous on large scales, and since the CMB is highly isotropic!! around us, this notional 
operator must reflect these observations by acting on the general metric of Eqn. ([T]) to yield a 
FLRW metric. Crucially, the scale factor of this FLRW metric must (within our assumptions) 
be given by Buchert's functional av{t), where the domain size to be used will be specified 
shortly. The argument leading to this proposal is as follows. 

We begin by making precise the notion of "homogeneity on large scales" , by requiring that 
in the range R = (Li,L2) of CGC sizes allowed by the fluid approximation validity criterion, 

*We proceed with the understanding that we will only be dealing with scalar quantities in a chosen time-slicing. 

^The CMB is isotropic after subtracting a dipole contribution due to the local motion of our galaxy. This 
motion can be thought of as a measure of inhomogeneities on small scales. The isotropy that we impose is 
therefore effectively an isotropy on large scales. 
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there exist a subset Rhom =C R satisfying the fohowing condition: When the microscopic 
matter density p{x,t) is averaged over a length scale L G Rhom, the resulting smoothed density 
must lose all information about the particular CGC over which the smoothing was performed. 
In other words, for L £ Rhom, at any given time t, the smoothed matter density must have the 
same value in all CGCs. It is in this sense that the smoothed matter density is homogeneous 
(see however, the discussion at the end of this section concerning an explicit smoothing of the 
metric). This definition can be related to the concept of statistical homogeneity if we assume 
that the total volume of the observable Universe is large enough that each CGC effectively 
samples one of an ensemble of configurations of matter distribution (see Ref. [26J for a careful 
discussion of this point). It is possible that homogeneity sets in at a scale too large to be 
considered as the size of the CGCs, namely that the range Rhom is not contained in the range 
R and that L'2 <^ -^^Hubbie is not satisfied. We will assume that this does not happen (for a 
more detailed discussion on this assumption, see Sec. 2.1 of Ref. [2]). We also make the crucial 
assumption that Rhom does not change with time. This will later allow us to consistently analyze 
light propagation using the effective FLRW metric. This assumption needs some justification 
and we will return to this issue in the discussion. 

Now define a coarse-graining cell to be an averaging domain D as used in Buchert's for- 
malism, with a characteristic comoving size L G Rhom- The averaged density is then {p)D(t), 
independent of the location of the particular CGC being considered. It also follows that the 
averaged expansion scalar {@)t> = 3{ax>/ax>) = (Vp/Vp) is homogeneous in the above sense 
(using Eqn. (|7ap) and describes the evolution of the physical volume of any given CGC. Thus, 
{@)v plays the role of the expansion scalar for a geometry whose infinitesimal spatial volume 
elements probe length scales of the size of our CGCs (see however the discussion at the end 
of this section). It is therefore natural to demand that (G)x) be the expansion scalar for our 
notionally smoothed out macroscopic metric. This is equivalent to demanding that the scale 
factor of this homogeneous and isotropic effective metric be identified with the functional aT>{t) 
(up to a constant factor which we choose to be unity by an appropriate choice of units for the 
effective comoving coordinates). We are thus left with an effective metric whose line element 
can be written as 

dsls = -dt^ + av{tf ( Y^A^ ^ r'^dn'^\ , (10) 

which probes comoving length scales of size L as chosen above. The 3-Ricci scalar of this 
effective metric is given by T^eff = Qkv/a^ and should be thought of as arising from the coarse 
graining of the geometry. Following Ellis' ideas from Sec. [2]then, 7^efr(= S Einsteini^^^TV)) is in 
general different from the averaged 3-Ricci scalar {TV)x> which arises from applying the matter 
averager to the inhomogeneous scalar curvature ((^)d = Smatti^^^T^))- The constant kx> has 
the usual interpretation as a parameter characterizing all possible homogeneous and isotropic 
3-metrics. We can write {TZ)v in terms of TZ^s by defining ATZ = {TZ)x> — 6kx>/a'jy. In doing so 
we are merely shifting focus from (the a priori unknown) {TZ)x) to another unknown ATZ. As 
the notation suggests, we think of ATZ as being a correction (a part of Pmac in the notation 
of Sec. [21 see Eqns. (jlip below) to the 3-Ricci scalar of the FLRW metric, although this 
correction need not be small compared to T^eff and need not evolve proportional to a^^. We 
choose not to normalize the magnitude of the constant kj), retaining instead the normalization 
of the functional a-D^t) at the initial time tj. The evolution equations appropriate at this 
scale are Buchert's equations ()7bp and (fTcj) together with the integrability condition ([9]). For 
completeness, and in the spirit of Ellis' arguments we can rewrite equations (jTbll and (iTc]) as 
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Gm,n.r. ^mac Pmac 



d-T) 4ttG , , 1 ^ , . 

— = ^(P)7? + i,Qv (lib) 

ax> o 6 



where we have symbolically denoted the specific components of the general modified equations 
to which the various terms of Buchert's equations belong. The integrability condition ([9]) which 
supplements these equations is purely a consequence of having non-zero corrections Pmac to the 
standard FLRW equations. As a quick check, we note that requiring the corrections Pmac to 
vanish implies both = (from Eqn. (jllbp ) and ATZ = 0, i.e. {TZ)v = 6A;x)/a|, (from Eqn. 
(jllap ): the integrability condition ([9]) is identically satisfied and we recover the usual FLRW 
solution. 

In the arguments leading to Eqns. (jlll) . we have ignored the effects of a virtual change 
in the averaging length scale L. Considerable effort has been spent by Buchert and Carfora 
|27j in studying such effects. In particular, in the spirit of the real space Renormalization 
Group formalism, these authors derive a novel curvature backreaction effect analogous to the 
kinematical backreaction Qx), together with a volume scaling effect. The inclusion of these 
effects would change the evolution equations satisfied by the scale factor a^i- We also note 
that Buchert and Carfora's techniques (involving Ricci flow of the inhomogeneous geometry 
into a smooth geometry) form one candidate for an explicit realization of the metric coarse 
graining operator Smet- Seen in this light, it is apparent that by accounting for the above 
effects, we would also have to appropriately modify our notion of large scale homogeneity (to 
include the volume scaling effect and the residual information from the curvature fluctuations). 
In the arguments leading to the effective metric (fTO]l . the scale factor ax> would be replaced 
by an appropriately "dressed" quantity characterizing the evolution of effectively infinitesimal 
volumes in the smoothed geometry. However, the main arguments leading to an effective metric 
which is FLRW, will not be altered given the existence of the range Rhom C R. For simplicity, 
we ignore these additional effects in this paper, and hope to investigate them in future work. 

It should be noted that a notion of coarse-graining is always implicit in standard FLRW 
cosmology. It is assumed that a coarse-graining (while satisfying the conditions of the fluid 
approximation) is possible on a length scale which coincides with the scale at which homogeneity 
is achieved (in other words, the existence of Rhom C R is assumed). Standard cosmology 
however ignores the possible modifications to FLRW cosmology brought about by the process 
of explicitly smoothing over the inhomogeneities - an issue which Buchert's approach addresses. 
Our key proposal in this paper is that the size of the averaging domain should be taken to 
coincide with the scale of homogeneity, and that (up to some corrections) the quantity ax>{t) 
corresponds to the scale factor of a template FLRW metric. 

We now illustrate how our proposal can facilitate the comparison of the predictions of the 
spatially averaged cosmology with observations, for example by computing the Di^{z) relation 
using the scale factor a^i (t) , which has a modified time evolution compared to the corresponding 
scale factor in standard cosmology. 
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5 Luminosity Distance using the Effective Metric 

In our approach, we take the effective or template metric to be FLRW, and since we are holding 
the size L of the averaging domain fixed, there is nothing conceptually new to be done while 
studying its null geodesies. The important difference, of course, is that the evolution of the 
scale factor is now given by the modified Friedmann equations (llip (subject to the caveats 
listed at the end of the previous section). We note that the use of the FLRW metric is, strictly 
speaking "unrealistic" in the sense that in the real Universe, light travels mostly in vacuum 
which has a non-zero Weyl tensor whereas in the FLRW cosmologies which are conformally 
flat, light travels through a geometry with vanishing Weyl tensor (see also Refs. [9] and [2j). 
This issue is clearly relevant in the standard approach as well. While we have not addressed 
this issue in the context of averaging, we emphasize that our approach must be thought of as a 
first step in gauging the effects of inhomogeneities, by using an explicitly averaged construction 
of the template FLRW model to be used. The arguments relating the redshift of a source to 
the scale factor ax) and to its comoving distance r from the observer, are exactly the same as 
in the standard approach. In particular, we have for the redshift z of a source that emits light 
at time t which is received by an observer at tim^ to j 

l + z = ^ (12) 



and the comoving distance r to the source is determined by solving the equation 

dr Z"*" dt' 



V^l - kvr^ Jt av{t' 



(13) 



A caveat to the above relations is that the "step-size" used in computing the integrals involved 
must be understood to be no smaller than the size L of the CGCs. Further, the entire con- 
struction must be interpreted as a "fitting template" that allows us to compare observations 
with theoretical predictions in a straight forward manner. With our choice of conventions, the 
magnitude of kf is a dynamically determined parameter of the model and determines the form 
of the luminosity distance Di^{z) as follows. Solving Eqn. for r = rf,m{z) gives 

rem{z) = Sk^{a{z)) ; a{z) ^ ^""^ 



Di^{z) = 075,(1 + z)rem(z) , (14) 
where Hd = dxi/ax) is the Hubble parameter and the function Skj, is defined as 

{{l/y/k^) sin (Vk^x) , /c© > 

X , kv = Q (15) 

(1/ \/—kx)) sinh {^/—kjyx) , < . 

In order to complete the picture we need the functional dependence of the Hubble parameter 
Hxi on the scale factor a© . This requirement is complicated by the fact that the system of 
equations (j7a|) . ()7bp . ([7cj) and ([9]) is only consistent, it does not close. In order to obtain the 



The subscript will indicate the value of the quantity at the present epoch to, and the subscript i, the value 
at the initial epoch ti. 
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required relation then, it is necessary to make certain assumptions about the evolution of the 
kinematical backreaction Qx> and the averaged 3-Ricci scalar {TZ)t>. Buchert et al. [21] have 
analyzed a class of scaling solutions of the form 

Qv = Qv,al ■ {R)v = nv,al. (16) 

The case n = —6, p = —2 is the only scaling solution with n ^ p. For the cases n = p, the 
backreaction is proportional to the averaged 3-curvatur^, 

Qv.=rnv. ; r = -^ ; n = -2'-±^ . (17) 

n + b 1 + r 



The forms of the integrability condition ([9]) and the solutions (jl6p indicate that the scaling 
solutions can be superposed to obtain new solutions. We will use two such superpositions to 
construct models which contain an accelerating scale factor aj) and which also yield analytic 
expressions for the luminosity distance Di^{z). The solutions we consider are of the form 

{n)T, = ^ + -al ; Qv = Pal, (18) 

in keeping with our earlier discussion of the quantity ATZ. Namely, we model ATZ by a single 
scaling solution. Here r and n are related as in Eqn. (|17p . and kv and /3 are constants with 
dimensions of (length) We are interested in models which admit an accelerating scale factor 
av, and will hence assume /? > and kv < (the reason for which will be apparent shortly). 
Inserting the relations (jlSp in Eqns. (jllap and (jllbp we find 

a^ 3 3 n + 2 

Hv 47rG Pi B ^ 

— = o 3" + q «o ' (19b) 

av 6 a^-, 6 

where we have defined pi = {p)v{tin) and hence written {p)v = Pi/a\y using the continuity 
equation (j7ap . The next two examples which we will consider are pathological, in that the 
luminosity distances Di^{z) constructed for these models are not defined for all z > 0. Never- 
theless, since these models yield analytical results, we shall display the expressions for Di^{z) 
that we obtain. In general the computation of Di^{z) expressions would have to be performed 
numerically. 



Case 1 : n = —3. The backreaction Qv in this model decays at the same rate as the averaged 
matter density {p)v- We choose initial conditions such that (3 > iirGpi, which ensures that the 
model accelerates (cf. Eqn. (|19bP ). This is consistent with Eqn. (|19aP due to the presence of 
the term kva'^ < 0, which allows the right hand side to be negative. Defining A = (2/3)(/? — 
A-nGpi) > we find 

a^Hv = avav = \/av [{-kv)av - A]^^'^ . (20) 



^ We have retained the notation of Ref. Ulj . The proportionality constant r should not be confused with the 
comoving coordinate in the effective metric. 
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Defining the "volume deceleration" qx> = —aT>/{ax>H^) and evaluating Eqns. (jl9ap and ()19bp 
at the present epoch to, we find 

_ ^mO _ Qvq _ 1 A ^ ^ _ 

kv = {avoHvo? (^mo - 1 - = - {avoHvo? (1 - 2qvo) < , (21) 

where Q^o = P/'^t>o present value of the backreaction, and we have defined the present 

matter density parameter Q^nO = (SttG/Oq)/ (3-ff|)g) = (SvrGyOj) / (3a|,^ff|,^) in the usual way. The 
function a{z) defined in Eqn. (jl4p which determines the luminosity distance -Dl(-z) reduces to 



a{z) = ^== In 
v—i^v 



(1 + VI - '^qvo) VTT^ 



Since gx)„ < 0, the redshift is constrained to take values < z < l/(— 2gx)o)- 



(22) 



Case 2 : n = —5/2. Here the backreaction Qi) decays slower than the matter density, and if 
we assume < f3 < AnGpi, then we have a situation wherein the effective scale factor initially 
decelerates, and then starts accelerating after a certain epoch which is determined by the values 
of the parameters /? and pi. Using Eqn. (jl9ap we have 

/SvrG 4/3 

Now we have = {Qmof^) - Qvo/i^H^^-,) as before with Qx^o = P/a^^, and 

- = (avoHvof (1 + ^^mo - 4gOo) ■ (24) 

The volume deceleration gx>o will be negative provided a-po corresponds to a sufficiently late 
epoch, namely if ax>o > {AirGpi/SP)'^ as can be seen from Eqn. ()19bp . and provided kxi < as 
seen from Eqn. (j24p . Assuming these conditions are met, we find for a{z) 



a{z) 



V-kv 
-In 



C^ - (0„o - 2qv,) VlT~z + C7 {1 - 2 (n„o - 2qv,) {VlT~z - l) + fi^o^}'^'] | 



C=^/l + nmO- ^qvo ■ (25) 

The range of allowed values for z is constrained in this case as well. For example, the quantity 
within the smaller braces in the previous equation must be positive, and this can be shown to 
imply that 

(z - z+)(z - z_) > 0, (26) 
where z± are the roots of the quadratic polynomial 

V{z) = Miz'^ + {2MiM2-l)z + Mf-l ; Mi = l + — ; - ^""^ 



2{^mO — 2qvo) 2 {Q,fno — 2qvo) 

(27) 
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An analysis of this condition reveals that unless 0,mO > ^Qvo^ range < Z- < z < must 
be excluded from consideration when qxiQ < 0. Requiring the arguments of the logarithms in 
Eqn. (j25p to be positive would also impose some condition on the allowed values for z. In any 
case, since there is no a priori reason to expect that > ^Qt>o holds, we see that this model 
also contains a pathology in general. 

Although both models considered above must be regarded as unphysical (since the luminos- 
ity distance of a source is a measurable quantity), consistent models also exist. For example, 
the case n = 0, r = — 1/3 corresponding to a constant backreaction, mimics the cosmological 
constant of the standard cosmology. Further, it may be argued that current observations seem 
to indicate that we live in a universe with kx) = and hence the above examples are not obser- 
vationally relevant in any case. To show that the formalism does admit consistent (accelerating) 
models with fep = as well, we now give an example of such a model. 



Case 3 : n = —1, kx> = 0. In this case, although the effective spatial curvature T^eff is zero, 
the average of the physical spatial curvature {TZ)v does not vanish, and is given by (see Eqn. 
(jlSp ) {TZ)t) = — 5/3/ax>- If we take /? > 0, which guarantees that the scale factor will accelerate 
for av > (iirGpi/ py^'^ (see Eqn. (|19bl) ). then we have < 0. The quantity a{z) is now 

given by 

r^o day SttG 2(3 

a(^) = / r-T72 ; = -^Pi . C2 = —. (28) 

Jar,{z) [Ciav + C2a%) ' 3 3 

Since Ci and C2 are positive, this integral is well behaved for all z. Further, since kx> = 0, 
we have Di^[z) = axi^iX + z)oi[z). Defining VLqq = {2/'i){QvJ H'^^) and noting that Ci = 
Q.rnoH'^^aj)^ and C2 = {2/^)QvQavo = ^QoH^^avo, we get 

, , (1 + z) dz 

Dl{z) = ^—^ / 77^ . 29 

Hvo Jo [l^„o(i + S)3 + J7go(i + 5)]V2 

In this notation, which is reminiscent of that in the standard LCDM model, if we use Eqn. 
()19ap evaluated at present time we also see that ^ImO + ^QO = Ij and that acceleration begins 
at a redshift given by 

As an example, in fig. [2] we have plotted the behaviour of D\^{z) in units of cH^^ in this 
model (for redshifts between z = and z = 1), setting $7^0 = 0.3, Qqq = 0.7. In this case 
the acceleration redshift is Zacc — 0.53. For comparison we also show the Di^{z) curve in the 
standard LCDM model with = 0.3 and = 0.7, for which the acceleration redshift is 
Za.cc = i2^A/^m)^^^ — 1 ~ 0.67. [Interestingly, a model in which the backreaction Qd behaves 
as at early times, has recently been studied by Li and Schwarz [28] in the context of spatial 
averaging of a perturbed FLRW Cosmological model.] 



We have therefore demonstrated the existence of models of the backreaction which can 
admit an accelerating scale factor. It is important to keep in mind though, that the models 
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With backreaction Q_D 
Lambda-CDM 




Figure 2: The luminosity distance (in units of cHj^^^) versus redshift for the backreaction model 



with n 
with Qr 



-1 defined in Eqns. ([28]) and ([21 
0.3 and = 0.7 (dotted fine). 



(solid line) and for the standard LCDM model 



considered above were not physically well-motivated, they served only to demonstrate the con- 
struction of luminosity distances using the modified dynamics. A more realistic scenario would 
probably contain a backreaction which becomes significant only at late times, when structure 
formation has advanced sufficiently far. It will be an interesting exercise to numerically con- 
struct luminosity distances for such models and compare them with those of the concordance 
model. 



6 Discussion 

We have argued that for the purposes of Cosmological study, an explicit coarse graining of the 
matter density in the Universe must be performed using the formalism of spatial averaging 
developed by Buchert. Due to the observed large scale homogeneity of the matter density and 
the isotropy of the CMB radiation, the notionally smoothed out template metric obtained after 
this coarse graining must be described by a FLRW line element, with the scale factor being given 
by Buchert's functional 075 (t). The procedure for constructing luminosity distances is therefore 
identical to that in the standard approach, except that the evolution of the effective scale factor 
must be given by the modified FLRW equations ([7]) supplemented by the integrability condition 
dl]). (The scale factor and its evolution will be further modified by including perfect fiuid 
sources with non-zero pressure [SJ and the curvature backreaction and volume effects j27].) Such 
a construction of luminosity distances must be interpreted as a convenient way of comparing 
observations with theoretical models, while retaining the essence of the physical inhomogeneous 
geometry. It will be an interesting project to compare such luminosity distances in realistic 
backreaction scenarios with current Supernovae observations in an attempt to constrain the 
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magnitude and evolution of the backreaction. 

We end by discussing some additional issues concerning an explicit smoothing of the metric 
which we avoided in the present treatment, and also some potential problems with our con- 
struction. Firstly, we return to the issue of holding the comoving range Rhom, and in particular 
the smoothing scale L, constant. This is justified only so long as the condition L <C -Z^HubbicC-^) 
is satisfied, where -Z^Hubbiel-^) is the appropriate comoving scale associated with the observable 
Universe at redshift z. For example, if we choose a scale L which satisfies this condition at 
the present epoch, the same scale will not in general be valid at the last scattering epoch of 
the CMB. However, it seems reasonable that for a small enough redshift range, the condition 
L <^ -Z>Hubbio(-2) can be satisfied with a single constant L. For epochs sufficiently far back in the 
past, the scale L would have to be reduced, and the construction of the Di^{z) relation would 
no longer be valid. 

Further, a change in the size of the averaging domain would in general affect both the scale 
factor ax> and the "constant" kxi in a non-trivial way. In particular, any change in kx> would 
render the effective FLRW metric (llOp meaningless. Also, a change in kj) need not necessarily be 
caused by a change in the averaging domain size. It seems likely [29] that an explicit smoothing 
of the geometry using Buchert and Carfora's technique will yield a 3-space of different spatially 
constant curvatures at different times, this difference being independent of any change in the 
averaging domain size with time and being possibly incommensurate with our assumption of 
T^eff oc a'^. This could be interpreted in terms of a change in kj) with time, and in such a case, 
one would be dealing with a different FLRW template at different times. The effects of such a 
construction are not clear at present. Clearly this issue needs more careful study, and we hope 
to address this in the future. Of course, if kxi is made a function of time, care must be taken 
that the matter content is now so chosen that Einstein's equations continue to be satisfied, as 
pointed out by Rasanen [31] in response to [30]. 

Also, as pointed out by Rasanen [2], the use of a homogeneous and isotropic FLRW template 
will necessarily ignore physical effects such as a non-trivial shear in the underlying geometry. 
A potential problem with our template construction is that effects of the non-trivial scalar cur- 
vature (3)7^(f, t) on light propagation are also simplified to a proportionality to . Although 
the full scalar curvature will show up in the evolution of ax)(t) via its spatial average {TZ)x>, 
it is possible that this construction may inaccurately model the real Universe. A study of the 
evolution of the effective spatial curvature as mentioned in the preceding paragraph may help 
resolve this issue. 

There is also the related issue of degeneracy, namely that several underlying inhomoge- 
neous models may reproduce the same backreaction effects when averaged. For example, it was 
shown [T7| in the context of the LTB solutions, that models with a vanishing spatial curvature 
(the so called marginally bound models) in the inhomogeneous geometry also have a vanishing 
backreaction after averaging. One is now faced with the situation wherein interpreting say 
luminosity distance observations within the inhomogeneous geometry may distinguish two in- 
homogeneous models, but interpreting the data after averaging would not make this distinction 
and would hence perhaps rule out both models. This issue has been recently highlighted by En- 
qvist and Mattsson (the last paper in Ref. [TT]) who show that a marginally bound LTB model 
can fit the data from Type la SNe without any spatial averaging. Since the standard FLRW 
matter dominated model is also (a special, homogeneous case of) a marginal LTB model which 
doesn't fit the SNe data [13], the formalism after spatial averaging would rule out both models 
by demanding a non-trivial backreaction instead. However, since the LTB (and in fact any in- 
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homogeneous) models themselves face a similar degeneracy issue (for example several different 
LTB models can currently claim to fit the SNe data [H]), this matter cannot be considered as 
settled, and also deserves further attention. 

Lastly, we note that Buchert's averaging scheme itself has been criticized in the literature 
(see, e.g., Ref. [20]) on the grounds that the effects of averaging in a noncovariant manner as 
done in this scheme, are likely to be gauge artifacts. It should be noted that within second- 
order perturbation theory, gauge invariance of the averaged equations has been demonstrated 
in [16] and [28]. It is also possible to approach the problem of determining modifications to the 
standard Cosmological equations via a covariant method such as Zalaletdinov's Macroscopic 
Gravity [32] for example. Such an attempt has in fact been made recently and has led to 
effective Cosmological equations very similar to those of Buchert (see Ref. [33J). The proposal 
in the present paper however, is attractive to us due to its simplicity in implementation, and 
may be treated as a phenomenological model of the more rigorous equations derived in Ref. 
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